A SIMPLE APPROACH TO GLOBAL REGIME OF THE RANDOM 

MATRIX THEORY 



LEONID PASTUR 

Abstract. We discuss a method of the asymptotic computation of moments of 
the normalized eigenvalue counting measure of random matrices of large order. The 
method is based on the resolvent identity and on some formulas relating expecta- 
tions of certain matrix functions and the expectations including their derivatives or, 
equivalently, on some simple formulas of the perturbation theory. In the framework 
of this unique approach we obtain functional equations for the Stieltjes transforms 
of the limiting normalized eigenvalue counting measure and the bounds for the rate 
of convergence for the majority known random matrix ensembles. 



I. Introduction 

Random matrix theory is actively developing. Among numerous topics of the theory 
and its various applications those related to the asymptotic eigenvalue distribution 
of random matrices of large order are of considerable interest. An important role 
in this branch of the theory plays the eigenvalue counting measure defined for any 
Hermitian or real symmetric matrix 



j,k=i 



(1-1) M n = [M$ 

as follows 

(1.2) N n (A) = ^{\ ( r ) e A} 

where A is a Borel set of the real axis R and {Aj n) }™ = 1 are eigenvalues of M n . 

One distinguishes several large-n asymptotic regimes for the probability proper- 
ties of eigenvalues (see e.g. [3J). In this paper we deal with the global regime, defined 
by the requirement that the expectation E(iV Tl (A)) has well defined (i.e. not zero 
and not infinite) weak limit 

(1.3) N(A) = lim E(iV n (A)). 

n— >oo 

This limit is called the Integrated Density of States (IDS). We shall see below that 
explicit conditions to be in the global regime may look differently in different cases. 

The IDS is a quantity to be found and analyzed first in any random matrix 
study, because it enters in practically any problem and result of the theory. A new 
wave of interest to the global regime is motivated by recent studies of the free group 
factors of operator algebras known now as free probability theory H. 
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In this paper we present a simple approach of the study of the random measure 
( |1.2j ) in the global regime. The approach allows us to find the limit ( |1.3| ) in many 
interesting cases, to show that the sequence of random measures (|1.2|) converges to 
this nonrandom limit either in probability or even with probability 1 and to find 
bounds on the rate of the convergence. 

The paper is organized as follows. In Section 2 we consider most studied random 
matrix ensembles, in particular, the Gaussian and the Circular Ensembles. This 
allows us to explain the method by using the simple and known setting of these 
ensembles. We call these ensembles classical because, first, they are indeed classic 
objects of the theory, and second, because of the role of the classical orthogonal 
polynomials in their studies (although, we almost do not use this technique in the 
paper). Section 3 is devoted to ensembles whose probability distribution is invariant 
with respect to unitary or orthogonal transformations and whose study is motivated 
by the Quantum Field Theory. In Section 4 we present new results on the form of the 
limiting normalized eigenvalue counting measure of the sum of two Hermitian or real 
symmetric matrices randomly rotated one with respect to another. In Section 5 we 
first discuss ensembles with independent but not necessary Gaussian entries. These 
ensembles are known as the Wigner Ensembles. Then we consider the ensembles 
that can be represented as the sum of the rank one independent operators. This 
form generalizes the sample covariance matrices widely used in multivariate analysis 
and is also motivated by statistical mechanics. 

Most of results, presented in the paper are known, sometimes for decades. How- 
ever they were obtained by different and often rather complicated methods while in 
this paper we derive them in the framework of an unique approach, that we present 
in three slightly different versions, according to the case considered. We do not give 
here complete proofs of all presented results, but only outline basic moments of their 
proofs. The complete versions of the proofs will be published in |3J. 



2. Classical Ensembles 

2.1. Gaussian Ensembles. We start from the well known ensembles among which 
the Gaussian Ensembles (GE) are most known. We restrict ourself by the technically 
simplest of Gaussian Ensembles, consisting of Hermitian matrices and known as the 
Gaussian Unitary Ensemble (GUE) because its probability distribution 

(2.1) Pn(dM) = Z- 1 exp (-^ Tr M 2 ) dM 
is unitary invariant. Here 

(2.2) dM = Y[ dMjj Yl dReM jk dlmM jk 

is the "Lebesgue" measure on the space of n x n Hermitian matrices. It is often 
convenient to write 



(2.3) 



M n = n- l ' 2 W n 
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where now W n = {Wjk}™ k=1 can be considered as the left upper corner of the semi- 
infinite Hermitian matrix W = {Wjk\°° k=v whose entries are complex Gaussian ran- 
dom variables defined by the relations 

(2.4) E(W jk ) = 0, E(W jk W^) = 2w 2 S jl S km . 

Thus the probability space in this case consists of these matrices and has as a prob- 
ability measure the infinite product of the Gaussian measures defined by relations 

(El). 

The relations (ETSj) and (|2.4j) define the global regime in this case. 



Theorem 2.1. For the Gaussian Unitary Ensemble defined above the sequence of 
eigenvalue counting measures converges with probability 1 to the nonrandom 

measure 

1 

47TW 2 J AnMv^y^ 



(2.5) iV sc (A) = — - / V8w 2 - \ 2 d\ 



i.e. N sc has the density 

(2.6) p sc (A) = (4tt«; 2 )- V8w 2 - A 2 

concentrated on the interval [—2\/2w,2^/2w]. The convergence N n to N sc has to be 
understood as the weak convergence of measures. 

The theorem dates back in fact to Wigner (see [|IJ). We give below the two proofs 
of the theorem to illustrate two methods that can be used in rather general situation 
of dependent and not necessary Gaussian entries. Both proofs as well as other proofs 
in this paper are based on the study the Stieltjes transforms of measures instead of 
measures themselves. In the random matrix theory the Stieltjes transform was used 
for the first time in paper WM and since then is proved to be a rather efficient tool 



of the study of the global regime. 

Recall that the Stieltjes transform f(z) of a non-negative measure m(dX), m(R) = 
1, is the function of the complex variable z defined for all non-real z by the integral 

(27) /« = /t^ Im ^°' 

Here and below we use integrals without indicated limits denote integrals over whole 
real axis. f(z) is obviously analytic for non-real z and satisfies the conditions 

(2.8) Imf ■ Imz > 0,lmz ^ 0, swpy\f(iy)\ = 1. 

y>i 

It can be shown || that any function f(z) defined and analytic for non-real z and 
satisfying conditions (|2.8|) is the Stieltjes transform of a unique nonnegative and 
normalized to 1 measure m(dX) and that for any continuous function (p(X) with a 
compact support 

(2.9) / ip(X)m(dX) = lim- / <p(\)1mf(\ + ie)dX 
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Besides the one-to-one correspondence between measures and their Stieltjes trans- 
forms is continuous if one will consider the weak convergence of measures and the 
convergence of the Stieltjes transforms that is uniform on all compacts in C\R. 

The use of the Stieltjes transform in this context is based on the spectral theorem 
expressing the Stieltjes transform 

(2.10) g n (z) = / 

of the normalized eigenvalue counting measure ( |1.2|) of a matrix M n via its resolvent 

(2.11) G(z) = (M n - z)~\ Imz^O 
by the formula 

(2.12) g n (z) = - Tr G(z) 

n 

Our proof has basic ingredient the following 

Proposition 2.1. The Stieltjes transform g n (z) of the eigenvalue counting measure 
(Uj^ for the Gaussian Unitary Ensemble defined by ( \2.1\ ) or by t \2.3j ) and has 
the asymptotic properties for Imz > yo = 5w: 

(2.13) lim E(g n (z)) = f sc (z), 

n— +oo 

(2-14) E(| 7 „(z)| 2 )< C 



w 2 n 2 



where 

(2.15) ln (z)=g n (z)-E(g n (z)), 
f sc (z) is the unique solution of the equation 

(2.16) 2w 2 f 2 + zf + 1 = 

verifying condition ( \2.8j ), and we denote here and below by C numerical constants 
that may be different in different formulas. 

To prove the proposition we need the following elementary facts: 

(i) for any two matrices A and B 

(2.17) (B - z)- 1 = (A - z)- 1 -(B- z)-\B - A) (A - z)~\ 

(the resolvent identity); 

(ii) if £ is a complex valued Gaussian random variable defined by 

E(C) = E(C 2 ) = 0, E(|C| 2 ) = 2^ 2 

and (f(z, z) is a differentiate function polynomially growing at infinity and 
having the same property of its derivatives, then 

(2.18) E(C>(C, 0) = E(|C| 2 )E = 2^ 2 E ; 
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(iii) for the resolvent G(z) = (A — z)^ 1 of any Hermitian or real symmetric matrix 
A we have 

(2.19) \\G(z)\\ < \lmz\~\ \G jk (z)\ < \lmz\~ 1 

where Gjk(z),j, k = l...n are the matrix elements of the resolvent. 

Now using ( |2.17| ) for the pair B = M, A = and applying ( |2.18| ) we obtain the 



system of identities for the moments 

(2.20) m p {z\, z p ) = E(g n (zi)...g n (z p )) 
of the random function g n with non-real arguments Z\, z p 

1 2w 2 

(2.21) mi(^i) = rn 2 (z 1 z 1 ), 

Z\ Zi 



1 2w 2 
m p (z 1 ,...,z p ) = mp_i(z 2 , ...,%) m p+1 (z u Zi t z 2t ...,z p ) + r p (z 1 ,...,z p ),p> 2. 

Z\ Z\ 

where 



2w 2 P 1 

r p (zi, ...,z p ) = 5Z E ( _ Tr ( G n(zi)Gl(z q ))g n (z 2 )...gn(zq-i)gn(z q+ i)...gn(z p j, 



n*zi * — ' n 

q=2 



Assume now that |Imz 9 | > y,q = 1,... for some y > 0. Then relations ( 2.1 2| ) and 
(|2.19| ) imply the bounds 



(2-22) 1^)1 < _L < I 

|lm2;| y 

1 2w 2 p 
\m„\ < — , rJ < 



PI — n ' P — 9 

n z y p+1 



Following statistical mechanics (see e.g. 0]) we can treat system of identities ( |2.21| ) as 



a linear equation in the Banach space B of complex valued sequences m of functions 
m = {m p (zi, Zp)}^ =l equipped with the norm 

(2.23) 1 1 m 1 1 = sup ?7 P sup \m p (zi, z p )\ 

P^ 1 |Im 2 ,|>y,g=l,..p 

for some 77 > 0. The equation has the form 

(2.24) m = Am + b + r 
where 

. w , 2u; 2 2w 2 

(Am) p (zi,...,^) = m p -i(z 2 ,-,Vi) "Vn(2i, Zi, z 2 , ^ P -i),p > 2, 



(i 
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and b = {—z x , 0, ...), r = {r p } p ^ =1 . It is easy to show that optimal with respect to 

7] bound for the norm of A is < 2^ wy~ l for r\ = \/2w . Choosing say y > 5w 
we have uniformly in n that the norm of A is strictly less than 1, the vectors m,b 
and r belong to B and 

(2.25) IH|<-2- 

Thus equation ( p.24|) is uniquely soluble uniformly in n. It is easy to check that this 
equation with r replaced by zero has the factorized solution 

q=l 

where / verifies ( |2.16| ). Thus, in view of ( |2.25| ) we have 



p „ 
(2.26) \m p {z 1 ,...,z p )-Y[f{z q )\ < 



q=l 



w p n 2 



uniformly in 

(2.27) \lmz\ > 5w. 

In particular we obtain that 

lim mi(z) = f(z) 

c 

|m2(z, ~z) — mi(z)mi(~z)\ < 



n 2 w 2 



Recalling definition (|2.20|) of m p we see that the first relation implies that f(z) satisfies 
(|2.8| ). The unique solubility of ( p.!6| ) in this class can be easily verified. The second 
relation is just another form of Q2.14j) . The proposition is proved. 



Remark. As was mentioned the technique of the proof is similar to the tech- 
nique of the correlation equations of the statistical mechanics (the Kirkwood-Salzburg 
equations, the Montroll-Mayer equations, etc. |J) combined with the mean field ap- 
proximation also widely used in the statistical mechanics. The reason to have here an 
analogue of the mean field approximation regime is again similar to that of statistical 
mechanics: the entries of the GUE matrices are all of the same order of magnitude 
(see (j2.3|) - (j2.4|)), like interactions in the Curie- Weiss model. Thus to obtain a non- 
trivial (not zero and not infinite) limit ( |1.3| ) (an analogue of an extensive quantity 
per unit volume in the statistical mechanics) we have to introduce the n-dependent 
normalizing factor n~ l l 2 in ( j2.3| ). This fixes the global regime scaling but also leads 
to vanishing of the statistical correlations and to the factorization of the moments 
(see ( [2.14 ) and ( 2.26|) ) and to a nonlinear self-consistent equation determining the 
first moment that can be regarded as an analogue of the Curie- Weiss equation for 
the magnetization. 
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To prove Theorem 2.1 denote by N sc the measure ( |1.2j ) corresponding to f sc (z) 
via (|2.9| ). In view of Proposition 2.1 and the Borel-Cantelli lemma we have for any 
fixed z with \lmz\ > 5w the convergence of g n (z) to f S c{z) with probability 1 for 
any fixed z. Since any analytic function is uniquely determined by its values on a 
countable set having at least one accumulation point we have the convergence of g n to 
/ with probability 1 on any compact of the domain \lmz\ > 5w. Now by continuity 
of the correspondence between measures and their Stieltjes transforms we obtain 
that the measure N n converges weakly to N sc with probability 1. Solving explicitly 
equation ( 2.1(j| ) in the class ( |2.£| ) we find that 

1 



(2.28) 



fsc(z) 



5«r 



where the radical is defined by the condition that it behaves as z as z —>■ oo. By 
using the inversion formula (|2.9j) we obtain (JL3.5[) . Theorem 2.1 is proved. 



The method was proposed in and subsequently used in |7|, |], ||, |T0| to study 
a variety of problem of random matrix theory and its applications. A certain disad- 
vantage of the method is that it is rather tedious. An important simplification of the 



method was proposed by A.Khorunzhy [11]. We describe now this simpler version by 



giving another proof of the previous proposition. 



form 
(2.29) 



Rewrite the first two equations of the system ( |2.21| ) for z\ = z, z 2 

1 2w 2 



z in the 



E(0n(z)) 



(2.30) 



E(b«(z) 



1 



— E(<fc(z)) 



2w< 



E(g 2 n (z)g n (z))+r 2 (z). 



z z 

Expressing —\jz in the first term of the r.h.s. of ( |2.30|) from ( |2.29|) we obtain 

2w 2 



(2.31) 



E(|7»WI' 



-E(g 2 n (z) ln (z))+r 2 (z) 



where r ) n {z) is defined in ( |2.15| ). By using this definition we can rewrite the expecta- 
tion in the first term in the r.h.s. of this relation as E(g n (z)\'j n (z)\ 2 )+'E(g n (z))E(\j n (z) | 2 ). 
Hence in view of (|2.19|) 

2 



(2.32) 



Wn(zhn(z) 



< 



\lmz 



■E (hJz) 



By using this inequality and (|2.22| ) we get from ( |2.3 1|) 



(2.33) 



:i- 



Aw 2 
\Imz\ 



)V(\ln(z) 



< 



2w 2 



|Imz| 4 n 2 



Thus, we obtain the bound ( |2.14| ) on the variance of g n (z) under condition ( p.27| ). 
By using this bound we replace E(g 2 (z)) in ( |2.29| ) by E 2 (g n (z)) with an error of the 
order 0(l/n 2 ) uniformly in the domain ( |2.27| ). Now by using standard compactness 
arguments we can prove that any subsequence of the sequence Ei(g n (z)) converges 
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uniformly on compacts of the domain (|2.27|) to a solution of equation ( j2.16| ) satisfying 
(|2.8|). Since this solution is unique we obtain other assertions of the Proposition 2.1. 



Remarks, (i). Analogous results are valid for two other widely used ensembles, 
the Gaussian Orthogonal Ensemble (GOE) consisting of real symmetric matrices 
distributed according the real analogue of (|2.1|), and for the Gaussian Symplectic 
Ensemble (GSE), consisting of self-dual Hermitian matrices having also the Gaussian 
probability distribution (see [[!]] for definitions and properties). This allows us to 
write that the density of the semicircle law of all three cases is concentrated on the 
interval [— 2y/]3w, 2y//3w] and has the form 

(2.34) p SC/3 (A) = 7 J— y /Af3w 2 - A 2 , |A| < 2y/pw. 

2pnw z 

for the GOE (p = 1), the GUE {(3 = 2), and the GSE {(3 = 4) cases respectively, 
(ii) We can consider an ensemble of the more general form 

(2.35) H n = + M n 

where M n is as before and Hn^ is a matrix whose eigenvalue counting measure 
has a weak limit Nq as n — > oo. We denote by fo(z) the Stieltjes transform of Nq. 
In this case we use the resolvent formula ( 2.17] ) for B = H n and A = Hn^ and a 



natural extension of the above arguments. We obtain an analogue of Theorem 2.1 
in which the Stieltjes transform f(z) of the limiting counting measure N is a unique 
solution of the functional equation 

(2.36) f(z) = f (z + 2w 2 f(z)) 



belonging to the class (|2.8| ) . The equation defines the the deformed semicircle law 
(see for its properties). 

2.2. Laguerre ensemble. The ensemble is defined as 

(2.37) M n = -A n A* n 

n 

where the n x n matrix A n has the probability distribution (cf. Q2.1|) ) 

(2.38) P(dA) = Z~ l exp (-^ Tr AA*^j dA 

n 

dA = J J dKeAjklmAjk. 

j,k=i 

In other words the entries Ajk,j, k = 1, n of A n are independent complex Gaussian 
random variables defined by 

E(A jk ) = E(A%)=0, E(|A, fe | 2 )=2a 2 . 

Note that the matrix A is not Hermitian. The name of the ensemble is recent (see e.g. 
|T2"|j ) and is related to the fact that in the orthogonal polynomial approach [|T| one uses 
in the case of this ensemble the Laguerre polynomials (recall that in the case of the 
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Gaussian Unitary Ensemble one uses the Hermite polynomials). The ensemble models 
a generic positively defined matrix. The real symmetric version of the ensemble in 
which A is nxm matrix with statistically independent Gaussian entries is well known 
since the 30 's in the multivariate analysis as the Wishart distribution and describes 



the sample covariance matrix of m random Gaussian n-dimensional vectors [JT3 



Theorem 2.2. Let the Laguerre ensemble of random matrices be defined as above. 
Then its eigenvalue counting measure converges in probability 1 to the nonrandom 
measure of the form 



(2.39) N L (A) = —I \I^-Y^d\ 

The proof of the theorem follows the scheme of that of Theorem 2.1, that is it 
is based on an analogue of Proposition 2.1. In particular, we have here an analogue 
of the important inequality Q2.14j ). However, the analogue of ( |2.29| ) has the form 

E(g n (z)) = -- z -2a 2 E(gl). 

As a result, the corresponding quadratic equation is 2a 2 zf 2 + zf + 1 = (cf. (|2.16|) ). 
This leads to (|2~39l) . 

Remarks, (i). Analogous results are valid for the cases when the matrix A n is 
real or quaternion. Thus we obtain the general formula for the density of the limiting 
measure 



(2.40) ^ W = _L_^^ <A<4/V, 

that is concentrated on the interval [0,4/3a 2 ] for the orthogonal (J3 = 1), complex 
(/3 = 2) and quaternion (f3 = 4) cases respectively. 

(ii). We can consider an ensemble of more general form 

(2.41) H n = + M n 

where M n is as in (FE37D and H [ n ] is a matrix whose eigenvalue counting measure 
/Vn°' ) has a weak limit iVo as n —>■ oo. We denote by fo(z) the Stieltjes transform of 
iVo- In this case we use the resolvent formula fl2.17| ) for B = H n and A = Hn^ and a 
natural extension of the above arguments. We obtain an analogue of ( |2.36| ) in which 
the Stieltjes transform f(z) of the limiting counting measure iV is a unique solution 
of the functional equation 

2a 2 



(2-42) f(z) = /, 



1 + 2a 2 / 



belonging to the class (|2.8|). 

(iii) We can also consider a more general case when the random matrix A n is 
the nxm matrix with Gaussian i.i.d. entries. This case can be treated similarly. 
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We discuss this case in more detail in Section 5 considering arbitrary distributed 
independent entries. 

2.3. Circular Ensemble. The ensemble consists of n x n unitary matrices whose 
probability distribution is given by the normalized Haar measure on U(n). The 
ensemble was introduced by Dyson together with its orthogonal and symplectic ana- 
logues (see for references and results). We discuss below the simplest Circular 
Unitary Ensemble (CUE) but one can obtain similar results for two other ensembles. 

It is useful to write eigenvalues \j of the ensemble in the form \j = e l9j , < 
9j < 2n,j = 1, ...,n and to introduce the normalized counting measure (cf. ( |1.2|) ) 

(2.43) N n (A) = \{6 3 e A} 

where A is Borel set of [0, 2n). An analogue of the Stieltjes transform ( |2.7| ) for 
measures on the unit circle is the Herglotz transform Q 

(2.44) h(z) = J -^- m (d6), \z\ ± 1. 
Respective inversion formula is (cf. 



ip(9)m(d9) = lim — / (p(6)Reh(re ie )d6. 
r->i-o 2tt J 



We use here instead of ( |2.18| ) the differentiation formula 

(2.45) E((p'(M)AM) = E(<p'(M)MA) = 

valid for any C 1 function (p : U(n) — > C and any Hermitian matrix A. This for- 
mula and the spectral theorem for unitary matrices according to which the Herglotz 
transform h n (z) of the eigenvalue counting measure can be written as (cf. ([2.12|) ) 

(2.46) h n {z) = - Tr ^-±^ 

n U — z 

allow us to write the following relations for the moments of h n (z) : 

(2.47) E{h n {z)) = -1, \z\ ± 1, 



(2.48) E(|/i n (z)| 2 ) - \E{h n {z))\ 2 < — , \z\ < - 

n z 4 

The first relation shows that E(iV n (A)) = | A|/27r for all n. This is easy to understand 
because the Haar measure is shift invariant. The second relation plays the role of 
(|2.14| ). By using these relations and following the scheme of proof of Theorem 2.1 
one obtains 



Theorem 2.3. Consider the ensemble of unitary matrices distributed according to 
the Haar measure on U(n) (the CUE). Then the eigenvalue counting measure ( 2.4ty 
of the ensemble converge in probability to the uniform measure on the unit circle. 
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Remarks, (i). Analogous statements are valid also for the Circular Orthogonal 
Ensemble and for the Circular Symplectic Ensemble (see [[1|] for their definitions and 
properties). 

(ii) . Note that unlike Theorems 2.1 and 2.2 where we have the convergence 
with probability 1, Theorem 2.3 asserts only the convergence in probability, despite 
the bound ( [2.48 ). The reason is that in Theorems 2.1 and 2.2 we can consider W n 
and A n for all n as defined on the same probability space of realizations of the semi- 
infinite matrices W = {Wjk\Yk=\ an< ^ A = {Ajk\Yk=\ equipped with the infinite 
product Gaussian measure. It is clear that similar natural and simple embedding 
does not exist for unitary matrices. [] This case can be regarded as an analogue of 
the triangular array scheme of probability and the Theorem 2.3 is an analogue of the 
(Tchebyshev) law of large numbers, while the case of the Gaussian and the Laguerre 
ensembles is analogous to the scheme of infinite number of i.i.d. random variables and 
Theorems 2.1 and 2.2 are analogues of the strong law of large numbers. To deduce the 
convergence in probability of N n (A) for a fixed Borelien A from the convergence in 
probability its Herglotz transforms for a fixed z,\z\ < j one has to use the argument 
of Section 4 of |pfl . 

(iii) . The simple differentiation formula (|2.45 ) as well as its version ( f4.5|) below 
allows one to give a direct proof of the asymptotic freeness of unitary and diagonal 



matrices as n — > oo (see p[ for definitions and results and [16], [wj for some 
related recent results). Existing proofs are based on the representation of the Haar 
distributed unitary random matrices U as the phase in the polar decomposition of 
the Gaussian distributed random matrix X with complex i.i.d. entries and on the 
approximation of the phase by polynomials in X. Because of singularities of the 
polar decomposition representation U = X(XX*)~ l l 2 these proofs are not simple to 
implement in all details. The approach based on the formula ( |4.5| ), i.e. on the shift 
invariance of the Haar measure, seems more direct and simple (see |23|] and Remark 
(iv) of Section 4). 



3. Invariant Ensembles 

In this Section we discuss the random matrix ensembles defined by the proba- 
bility law 

(3.1) P{dM) = Z' 1 exp(-n Tr V{M))dM 

where M is a real symmetric or a Hermitian or a quaternion self-dual Hermitian 
matrix and V is a bounded below and growing sufficiently fast at infinity function of 
respective matrix. For V(X) = X 2 /Aw 2 we obtain the Gaussian Ensembles that were 
considered in the previous Section. In this paper we restrict ourselves by polynomial 
Vs. As it was in the case of the Gaussian Ensembles we discuss here the technically 
simplest Hermitian matrices. This subclass of ensembles ( |3.1| ) is motivated by Quan- 
tum Field Theory (see e.g. review \TB\). Following Quantum Field Theory we will 



1 Although one can always use the probability space that is the product over all n of the probability 
spaces consisting of the groups U (n) with the normalized Haar measure as the probability measure. 
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call V the potential. We will give below a result for convex Vs. More general case 
see in [|J. 

Theorem 3.1. Consider the random matrix ensemble consisting of Hermitian n xn 
matrices distributed according to ( \3. 1\) in which the potential V is a convex polynomial 
of an even degree 2p and growing at infinity. Then the eigenvalue counting measure 
of the ensemble converges in probability to the nonrandom measure whose density 
is concentrated on the interval (a, b) and has the form 

(3.2) p{\) =p 2p _ 2 (\) v /(b- A)(A-a), a<A<6 
where 

(3.3) MW 4fM * 

is a positive on (a, b) polynomial of the degree 2p — 2 and a and b are uniquely defined 
by the equations 

(3-4) / ^ \7 g = 2n5 lq , q = 0,l. 



The proof of the theorem follows again the scheme of the proof of Theorem 2.1. 
In particular we have the analogue of Proposition 2.1 

Proposition 3.1. Under the conditions of the preceding theorem the Stieltjes trans- 
form g n (z) of the eigenvalue counting measure of the ensemble ( \3. 1\ ) has the following 
properties for \Imz\ > y and a certain y depending on V 

(3.5) hm E(g(z)) = f(z), 

71— »00 

const 



(3-6) E(\gi(z)\) - \V(g n (z))\ < 



n 



2 



where f(z) is a unique solution of the quadratic equation 
(3.7) f + V '(z)f + Q(z) = 

satisfying ( \2.$) and 

(3.8, 



z — ii 

in which N(dfi) is the measure corresponding to f. 



To prove the proposition we use the differentiation formula |T5 

(3.9) E<y(M) • B) - nE(>(M) Tr V'(M)B) = 

valid for the matrix distribution ( |3.1| ), any function ip : R — > C whose derivative 
is polynomially bounded on the whole real line and any Hermitian matrix B. By 



2 See Remark (ii) after Theorem 2.3. 
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applying this formula to the matrix element of the resolvent G(z) = (M — z)^ 1 we 
obtain the relation 

(3.10) E(^(z)) + E(i Tr G{z)V\M)) = 

n 

where g n (z) is defined in (|2.12|) . The identity 

(3.11) G{z)V\M) = G{z)V\z) + G{z){V\M) - V\z)) 
allows us to rewrite this relation in the form 

(3.12) E(g 2 n (z)) + V'{z)B(g n (z)) + Q n (z) = 
where Q n {z) is defined as 

(3.13) Q n (z) = E(- Tr Q n (z, M)), Q n (z, M) = G{z){V'{M) - V'{z)) 

n 

and is a polynomial of the degree 2p — 2 if V(z) is a polynomial of the degree 2p. 
It is easy to see that for V(z) = z 2 /Aw 2 (|3~lp coincides with (^Wj) . Thus ( g3g) is 
an extension of ( |2.16| ) to the more general case of distribution ( |3.1| ), where V(X) is a 
polynomial of an even degree bigger than 2. An analogue of Q2.33D has the form 

(3.14) 

^ " Ti 17o) E (l^ z )D ^ rf7^ E (- Tr ^(z)q n (z, M)) + 1 

|lm2; • k (2;J| |k '(2;J| n n^|(lm2;) d l/ 

where 

q n (z, M) = Q n (z, M) - E(Q n (z, M)). 

In the Gaussian case the first term in the r.h.s. of ( |3.14j) is absent. 
In view of the inequality 

|E(- Tr ln (z)q n (z,M))\ < E(| 7n (z)| 2 ) 1 / 2 E(i Tr q n (z, M)q*(z, M)) 1 ^ 2 
n n 

where 7 n (z) is defined in (|2.15|) , it seems that the most natural way to obtain ( |3.6| ) 
is to prove the estimate 

,1 , . „ . N N const 

E(— Tr q n (z,M)q (z,M)) < — — 
n n z 

thereby reducing the estimation of the variance of g n (z) to that of q n (z, M). Unfor- 
tunately, we do not know the proof of the last estimate based on the differentiation 
formula ( |3.9|) and similar to that in the second proof of Proposition 2.1. Thus we 
refer the reader to works [l!|,[20|, where the bound (|3.6| ) is proven for all Imz| > 
by using a combination of the orthogonal polynomials and variational techniques. 
A simple proof of a weaker version of ( |3.6| ) with n instead n 2 in the r.h.s. will be 
given in also by using the orthogonal polynomial technique. Any of these bounds 
allows us to replace E(g n (z) 2 ) in ( |3~5| ) by E(g n (z) 2 ) = E(g n (z)) 2 = f 2 (z). Besides, 
by applying (|3.9[) to (f(M) = M, we obtain the equality 



(3.15) E(- Tr MV'(M)) 

n 
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that allows us to prove that all coefficients of the polynomial Q n (z) are uniformly 
bounded in n . After that simple compactness arguments yield that the limit of any 
convergent subsequence f nj { z ) satisfies Q3.7|) . 

Having Proposition 3.1 we can prove Theorem 3.1 by the following arguments. 
By solving the quadratic equation ( |3.7| ) we find that the measure N has the bounded 
Holder density p, that for a convex V the support of the measure N corresponding 
to / is a finite interval (a, b) and that (see ]3[]) 

[ b p(u)du V'(X) , , „ 

(3.16) v.p. / ^W^ = ^ AG (a, 6) 

where the symbol v.p. J denotes the singular Cauchy integral. Regarding this relation 
as a singular integral equation for p(A) and using standard facts of the theory of 
singular integral equations |n] we find that the bounded solution of the equation 
has the form 

(3.17) P (X) = \VW) f YMz^M ^_ R(X) = (b-X)(X-a), 

^ J a p. — A y/R(fJL) 

provided that 

f b V'(u)du 

(3.18) / , ! = = 0. 



y/Rdij 

This gives condition ( |3.4j ) for q = 0. Besides we have the normalization condition 



(3.19) / p(fj)dfi = 1 

J a 

that can be rewritten in the form (|3.4j) for q = 1 by using ( p.l7|) . It is clear that the 
integral is positive for a convex V and that it is a polynomial of degree 2p — 2, if V 
is a polynomial of the degree 2p. The unique solubility of system ( |3.7|) can be proved 
by using the implicit function theorem ||. 

Remark. Consider the case of the monomial V(A) = \X\ 2p /2p. In this case 
above formulas can be written in the form 

27tl 2 p-i J\x\ Vt - A 2 h P 

where 

t a dt 

1 n 



1 CX. 



These formulas are also valid for non-integer p, i.e. for potentials of the form V(A) = 
|A| a /a provided that a > 2. For this case the formulas were obtained in p2| by 
another method. They can also be obtained by a version of the method presented in 
this Section. In this version we use the identity 

(3.21) G{z)V'{M) =G{z)V'{X) +G{z){V'{M) -V'{X)), z = X + ie 

instead of fl3.11| ). It can be shown that this allows us to obtain the final formulas 
^)-( p.4[ ) for example for non-polynomial (and even non- analytic) Vs provided that 
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they are convex, even, grow faster than logarithmically at infinity and are of the class 
C 2 on any finite interval (—L,L). 



4. Law of Addition of Random Matrices 
Consider Hermitian matrices of the form 

(4.1) V n A n V* + U n B n u: 

where V n and U n are random independent unitary matrices distributed both according 
to the normalized Haar measure on U(n), A n and B n are Hermitian matrices such 
that their normalized eigenvalue counting measures N An and N Bn converge weakly 
to the limits N A and N B respectively and satisfy the condition 

(4.2) sup / \\\N An , Bn (d\) < oo. 

n J 



Theorem 4.1. The normalized eigenvalue counting measure of the ensemble of ran- 
dom matrices defined above tends in probability as n — > oo to the nonrandom limit 
whose Stieltjes transform f(z) is a unique solution of the system of functional equa- 
tions 

f(z) = f A (z-A B (z)f- 1 (z)) 

(4.3) f(z) = f B (z-A A (z)f-\z)) 

zf(z) = A A (z) + A B {z) - 1 

where f belongs to the class ( $.8{ ) and A A and A B are analytic for non-real z and 
such that 

(4.4) swpy\A A (iy)\ < oo, swpy\A B (iy)\ < oo 

y>i y>i 



The theorem was proved in |25[ for the case of uniformly bounded in n matrices 
A n and B n by computing asymptotic form of moments of the sum via the moments 
of summands. This requires rather involved combinatorial analysis and impose the 
boundedness condition on matrices. 

We outline now the proof J23| based on the same ideas as above, i.e. on the 



resolvent identity and on a certain differentiation formula. The formula used in this 
case is 

(4.5) E(<p'(UBU*)[UBU*,C\) = 

where (p : R — > C is a C 1 function whose derivative is polynomially bounded on the 
real line, B and C are Hermitian matrices, [B, C] = BC — CB and the symbol E(...) 
denotes the integration over U(n) with respect to the Haar measure normalized to 1. 
The formula can be easily derived from the shift invariance of the Haar measure. 

Assume first that the norms of the matrices A n and B n are bounded uniformly 
in n. By applying (O) to the resolvent identity ( j2.17| ) relating the resolvent G of 
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matrix ( |4.S| ) and the resolvent G\ of matrix A n we obtain the matrix identity 

E(G— Tr G) = GiE(- Tr G) - G{E(G— Tr GUBU*). 
n n n 

Assuming now asymptotic vanishing of the fluctuations of normalized traces, that we 

had in all cases above, we can rewrite this matrix identity in the form 

(4.6) E(C) = Gi ( z - A /"/^ + o(l), rwoo 



/n(z) 



where 



(4.7) A Bn = E(— Tr U n B n U*G), f n = E(— Tr G). 

n n 

and |Inx2r| is large enough to guarantee inversibility of the argument of Gi uniformly 
in n. Applying to (|4.6|) the operation Tr we obtain the prelimit form of the first 
equation of system ( |4.3| ). The second equation follows from the analogous procedure 
in which the roles of A n and B n are interchanged (recall that U n and U* have the same 
distribution). The third equation is the limiting form of the identity n -1 Tr G(z — 
A n — U n B n U*) = 1 and of ( |4.7| ). The unique solubility of ( [4.3| ) follows from the 
implicit function theorem applicable for large |Imz|. The proof of the vanishing of 
the correlations, more precisely, a bound similar to ( 2.14 ), is based on the same idea 



sec 



To obtain the general case ( |4.2| ) we truncate eigenvalues of A n and B n by a large 
number T and use the minimax principle to control this procedure as T — > oo, the 
compactness arguments and the unique solubility of (|4.3|) in the class 



Remarks, (i). Since the normalized eigenvalue counting measure is unitary 
invariant and the Haar measure is shift invariant we can restrict ourselves without 
loss of generality to matrices of the simpler form 

(4.8) H n = A n + U n B n U* n 

This form have, for example, the matrices (|2.35|) of the deformed GUE. Indeed, 



any matrix belonging to the GUE can be written in the form f n A n $*, where \l/ n 
is the matrix of its eigenvectors, distributed uniformly over the U(n) according to 
the Haar measure, A n is the random diagonal matrix of eigenvalues and \l/ n and A n 
are independent jQJ. Besides, according to Theorem 2.1 the normalized eigenvalue 



counting measure of A n converges with probability 1 to the semicircle law (2J3). Thus 
\I/ n plays the role of U n , A n plays the role of B n and N B is given by ( |2.5| ). It can be 
easily checked that in this case the system (|4.3|) reduces to ( |2.36| ). Analogous fact is 
also valid for the deformed Laguerre ensemble (|2.41 ) and also for certain classes of 
random operators acting in l 2 (Z d ) ||26|| . 



(ii). It can be shown that the theorem is also valid in the case when matrices 
A n and B n in ( |4.1| ) are also random, but independent of U n and V n and N^ n and Nb„ 
converge weakly in probability to the nonrandom Na and A?#- Then the ensemble 
of deformed covariance matrices ( |2.36| ) in which the random vectors a\ are uniformly 
distributed over the unit sphere in C n also has form (O). As for the form ( [4. 1[ ) , 



it is the case for the sum of two independent matrices distributed each according to 
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the law ( p.l| ) with possibly different polynomials V\p,- In this case the condition ( |4.2| ) 
follows from ( |3.15| ). This case was considered in |23| by using formal perturbation 



theory around the Gaussian ensemble. Thus, we see that Theorem 4.1 describes in 
a rather general setting the result of deformation of a random matrix by another 
matrix randomly rotated with respect to the first and allows us to find the limiting 
eigenvalue counting measure of the the sum (of the deformed or of the perturbed 
ensemble) provided that we know these measure for the both terms of the sum. 



(iii). For any function f(z) satisfying ( |2.8|) one can introduce the "selfenergy" 
E(z) by the relation 

(4.9) f(z) = -(z + E(z))-\ 

It can be shown that S(z) is also analytic for non-real z, and has the same property 
ImS(2;)Im2; > 0, Imz ^ as / (see ( p.8|) ). Denote by z(f) the functional inverse 
of f(z) and set £(z) = R(f)- Then it is easy to see that ( (4.3| ) is equivalent to the 
relation 

(4.10) R(f) = R A (f) + R B (f) 

where Ra(I) and -Rs(/) are the selfenergies corresponding to /a(-z) and /b(^). The 
inverses of the Stieltjes transforms of limiting eigenvalue counting measures were used 
14]] in the qualitative study of the support of of these measures in the case (|5.16|) 



m 



below where 

ta(dt) 

ITT/' 

Relation ( fTTcD was noted in |f28[ for the case when Na and Nb are both the semicircle 



(4.11) R(f) 



laws (|2.5|) , when R(f) = w 2 f. The general form of this relation was proposed by 
D.Voiculescu in the context of the operator- algebras theory and its new branch known 
as the free probability theory (see for results and references). In this theory the 
semicircle law plays the role of the Gaussian distribution and the measure defined 
by (|2.4CI|) (more generally, by formula ( |5.16| ) below) plays the role of the Poisson 
distribution. 

(iv). Similar technique can be applied to multiplicative families of positive 
defined Hermitian and or unitary matrices and gives results P7j that generalize and 



simplify those of |T|| and also gives a more direct proof of certain results obtained 



for these ensembles in the context of free probability . 



5. Matrices With Independent and Weakly Dependent Entries. Tiny 

Perturbations 

5.1. Wigner Ensemble. The proofs outlined in previous Sections for simplest ar- 
chetypal ensembles, the GUE first of all, can be elaborated and used in rather gen- 
eral case of Hermitian, real symmetric or self-dual Hermitian random matrices whose 
entries are independent or weakly dependent modulo symmetry conditions. The ma- 
trices can be written in the form (cf.( |2.3| )) 

(5.1) M n = n~ 1 / 2 W ri 
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where matrix elements w!£ of the matrix W n still satisfy (|2.4j) but their probability 



jk 

laws P^\dW) are not necessary Gaussian and may be n- dependent. We call these 
ensembles the Wigner Ensembles. In this general case we have to use instead of 
differentiation formula (|2.18|) the formula |31[ 



(5-2) E(^(0) = E ^ E (^ (a) (0) + e 



a! 

a=l 



where K a are semi-invariants (cumulants) of a real-valued random variable £, (p : 
R -> C is a function of the class C s+1 and |e,| < C s sup a . |E(|£| S+1 ). 

Another version of the method is based on the perturbation expansion of matrix 
elements of the resolvent in a particular matrix element of the matrix M n . Indeed, 
according to Section 2 an important moment of the method is the asymptotical 
computation of the expectation 

(5-3) E(i=f>rf) 

v k=l 

basing on various differentiation formulas (see formulas ( |2.18|) , (|2.45|) , (|3.9| ), and( [5.2|) 
above). However, in the case of independent entries satisfying Q2.4|) , this requires the 



knowledge of dependence of Gjk on up to linear terms only. Indeed, writing the 
resolvent identity (|2.17|) for B = rr x l 2 W n and A = n~ x l 2 W \ w (-n) =Q we obtain 



(5.4) G jk = G% - n-^(W^Gf k Gi + W^GJgS) + r n 

where 

\W in) \ 2 

(5-5) G kl = G\ w(n)=0 , |r n |<y^. 

vv ki ~ u nllmzr 



Substituting ( f>A\) in ( ^.3|) and taking into account that G kl is independent of we 
can perform explicitly the expectation with respect to and obtain the relation 

2 n 

(5.6) E( 9n ( Z )) = -l/,-^E£GX) + 0( max E(\W^\ 3 ) /n^llraz] 4 ). 

ZTl — ' 0<i,k<n J 

k=l ~ 

Now we can use ( |5.4|) in the opposite direction to replace matrix element of G kl by 
those of G. Thus, if 

(5.7) sup max EflW^I 3 ) < w 3 < oo, 

n 0<j,k<n ■> 

we obtain the analogue of (|2.29|) in the case of independent entries satisfying ( p.4|) 
and ( ^.7|) with the error of the order n -1 / 2 . Similar arguments allows one to prove an 
analogue of ( |2.14| ) with the r.h.s. of the order n~ 1//2 . This is sufficient for the proof 
of an analogue of Theorem 2.1 for the independent entries satisfying ( |2.4| ) and ( p.7| ) 
and with convergence in probability instead of convergence with probability 1 (see 
also Remark (ii) after Theorem 2.3). 
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We list now several recent results obtained by combinations of approaches based 
on formulas (|5.2| ) and ( |5T4| ) (for an account of previous results see p^). 

(i) Semicircle Law. The normalized eigenvalue counting measure converges 
weakly in probability to the semicircle law ( |2.5| ) if and only if in addition to ( |2.4| ) for 
any r > matrix elements (|5.1|) of satisfy the condition 



(5.8) 



\z\ a P$\dW) = 0, 



reminiscent the well known Lindeberg condition of the validity of the central limit 
theorem. This fact is known since the seventies, see [28|] for the sufficiency of some- 
what stronger version of ( |5.8|) and [3(J for the necessity and sufficiency of Q5.8|). 
However these results were obtained by rather complicated method. In || we give a 
simple proof based on the approach of this paper. 

(ii) 1/n expansion JJT]. By using Q5.2| ) and assuming that {W^} are identi- 
cally distributed (modulo symmetry conditions as usually) and have s + 1 finite mo- 
ments one can construct 1/n - expansion of moments ( [2.20|) in powers of l/n l , I < s, 
with the error of the order l/n s+1 / 2 provided that the complex spectral parame- 
ter z verifies ( ^.27| ). We give here the two results for the real symmetric matrices 



M£ = n-Va(l + ^)W», E(H/ ( 



k 



0, E{Wj k ) 



w 2 and 



supE(|^| 5 ) 



< oo. 



We have then the following asymptotic formulas: 



mAz) 



fs 



1 + 



n 



+ 



(i-w 2 fi(z)y i-w*fi(z) 



where a 



E(\W jk 



(")|4^ 



0{n 



is known as the excess of random variable W. 



l )-3E(|W^ 

and is assumed to be independent of n, and f sc is defined in (|2.28| ); 



(n) 
jk 



m 2 (zi, z 2 ) - m 1 (zi)m 1 (z 2 ) = n 2 c(zi, z 2 ) + 0(ri 



where 
(5.9) 



2w 2 



:i-^m Zl ))(i-w 2 fUz 2 )) 



w 



fsc(Zl) — fsc{z 2 ) 



Z\ - Z 2 



+ *£(*l)/£(*2) 



(iii) Central Limit Theorem Assume in addition to (O) that the forth 



moments of exist and are independent of j, k and n. Then for z from the domain 



jk 



Q2.27] ) the random function g n (z) — E(g n (z)) converges in distribution to the Gaussian 
random function with zero mean and the covariance (|5.9| ). 



5.2. Sample Covariance Matrices. In this Subsection we consider an ensemble of 
random matrices, generalizing the Laguerre ensemble of Subsection 2.2 and its real 
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symmetric version known as the Wishart Ensemble of the sample covariance matrices. 
||13|| . Respective matrices have the form 



(5.10) AI mn A mn T m A m n 

n 

where A TO>n are nxm random matrices whose entries A^' n ^ are i.i.d. complex random 
variables satisfying conditions 

(5.11) E(4T>) = E((4r } ) 2 ) = 0, E(\A^) = 1, 

(5.12) sup max E(\A ( ™' n) \ 4 ) < a 4 < oo, 

m,n l<j<n,l<k<m 

and T n is a diagonal matrix. We assume that 

Tfl 

(5.13) m — > oo, n — > oo, > c < oo, 

n 

and that the normalized counting measure 

(5.14) <r m (A) = h{h e A} 

n 

of eigenvalues ti, I = 1, m of T n has a weak limit 

(5.15) o" m (A) — > cr(A), m — > oo. 

In particular, t^,/ = 1, ...,m may be i.i.d. random variables independent of A m ^ n . 

Theorem 5.1. Under conditions listed above the eigenvalue counting measure of ma- 
trices ( \5.1Q ) converges weakly in probability to the nonrandom measure whose Stieltjes 
transform is a unique solution of the functional equation 

ta(dt) 
T+tflz) 



(5.16) f(z) = -lz-c 

in the class (WW. 

We outline the scheme of the proof, based on the same idea as above, i.e. on 
the careful analysis of the result of infinitesimal as n — > oo changes of respective 
matrices. 

Start again from the resolvent identity ( 2.1 7| ) written for the pair B = M n , 
A = 0. Applying to the identity the operation E(n~ 1 Tr ...) we obtain 



(5.17) E(g n (z)) = -\- -L £ tjE ( £ A^G jk A 



(m,n) \ 

zn* — ' ' * — ' / ' / 

1=1 j,k=l 



We can use now the scheme of proof of Theorem 2.1 using ( |5.2j ) instead ( |2.18| ). It is 
more convenient however to apply here a somewhat different scheme. It is analogous 
to that based on relations (|5.3|), (|5.4|), and ( |5 . 5| ) in the case of the Wigner Ensembles 
of the previous subsection, however applied not to individual matrix elements but to 



A SIMPLE APPROACH TO GLOBAL REGIME OF THE RANDOM MATRIX THEORY 21 

the columns a/ = {n~ l / 2 A^ }™ =1 , 1 = 1, m of the random matrix A^ m,n \ Treating 
the columns as vectors of C n we can rewrite ( |5.17| ) as follows 

(5.18) E(9n{z)) = ---—Y] *iE((Ga,,a,)) 

z zn ' 
«=i 

where (., .) is the scalar product in C n . Since vectors a\ are independent we perform 
first the asymptotic computation of the expectation with respect to a\ in the Z-th 
term of the sum like we did in the previous subsection for . To this end we use 
the formula, giving in the explicit form the result of perturbation of the resolvent Gc 
of an arbitrary matrix C by the rank one matrix L a , a G C n defined by its action on 
any vector x G C n as L a x = (x, a) a: 

(5.19) (C + L a - z)- 1 = G C - G c L a G c (l + (G c a, a))" 1 . 

The formula can be easily derived from the general resolvent identity fl2.17|) . By 
applying the formula to C = M n \ ai=Q we obtain that 

ti(Gia h ai) 



1 + ti(Gia h ai) 



(Ga h ai) - 

where Gi = G| ai= o. This relation will play the role of ( |5.4j ). Indeed, assume first that 
for some finite T and ct2 

(5.20) sup max < T, sup max ||a/|| < a2 

n l<l<m n l<l<m 

Since Gi does not depend on a; and since random vectors {a/}[^ 1 are mutually inde- 
pendent one can find from ( |5.11| ) and (|5.12|) that 

E l ((G l a l ,a l )) = -TrG l =g^ 
n 

|E,((G ia „ ai ))-iTYG,|'< -g^ 

where the symbol E/(...) denotes the operation of the expectation with respect the 
vector ai only. 

These relations allow us to present ( |5.18|) in the form (cf. 



^ m 

(5.21) E( Sn ) = y)tiE 

where now 



a [l) 

(Jn 



1 + 



W <£(i±M, |Imz| > ;/„, 
n\Lmz\ z 

and ?/o depends on T and on a 2 of (|5.20|) . Besides, applying to (|5.19|) the operation 
- Tr ... we obtain that 

n 

_ (i) _ ti(Gfai, a t ) 

n 1 + ti(Gia h ai) 
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\9n{z) 



1 



By using three last relations we can write instead (|5.21|) for |Imz| > yo and n— > oo 
(cf. O) 



(5.22) 



E(0„(z)) 



E 



tg n {z) 



<r m (dt) + o(l) 



1 + tg n (z) 

where a m (dt) is defined in ( |5.14p . This is an analogue of ( |2.29| ). Similar arguments 
allow us to prove an analogue of fl2.14| ). As a result we obtain ( |5.16| ) in the case 
of bounded ti and a/. General case can be obtained from the proven one by using 
the the analyticity of the Stieltjes transform up to the real axis, the truncation of ti 
and ai, the minimax principle to control the truncation procedure, the compactness 
arguments, and the unique solubility of (|5.16|) in the class Q2.8Q . The latter results 
from the the implicit function theorem. 

Remarks, (i). In the case when m = n and a(dt) has only one atom at t — 2a 2 
we obtain (|2]42|). 

(ii). Similar arguments shows that the eigenvalue counting measure of deformed 
ensemble ( 5.1Up 

(5.23) H n = Iff + M m , n 

where M m>n is defined by ( |5.10D and iff has the limiting eigenvalue counting measure 
iVo (like in (|2.35|) and in ( |2.41| ) also tends weakly in probability to the nonrandom 
limit whose Stieltjes transform is a unique solution of the functional equation 

ta(dt) 
l+tf(z) 

This functional equation was derived first in by another and rather complicated 
method. The method was based on the study of the sequence of matrices iff, p = 
1, .... m defined as 



(5.24) 



/(*) = fo\z-c 



(5.25) 



(m,n) ^ (m,n) 
Ik 



1=1 



and "interpolating" between iff and (Hn l> ) = H n . Asymptotic computation of the 



differences Tr (Hn 



(p+i) 



Tr 



(p) 



_1 based on formula (|5.19|) lead to the 



first order partial differential equation for 



f(t,z) 



1 



lim - Tr (iff -z)' 1 . 



Solving the differential equation subject the conditions /(0, z) = fo(z), f(l, z) = f(z) 
one obtains (|5.24 ). 

The derivations of equation ( |5.24j ) given later in 



Til and as well as the 



proof outlined above are more simple and direct. On the other hand, the sequence 
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H^ P) in (]p5D can be regarded as a matrix version of the sum of independent random 
variables with varying upper limit used often in the study of limit theorems and 
processes with independent increments. Similar observation was used recently in |2| 
to construct free (non-commutative) analogues of these processes where, in particular, 
an analogous partial differential equation was obtained (called in the complex 
Burgers equation). 



Acknowledgment 

I am thankful to A.Khorunzhy, B.Khoruzhenko, A.Stoyanovich and V.Vasilchuk 
with whom the results mentioned above were obtained. The final version of this paper 
was written when I was participating the semester "Random Matrices and their 
Applications" at the MSRI. I am grateful to Prof.D.Eisenbud for the hospitality at 
the Institute, and to organizers of the semester P.Bleher and A. Its for the invitation. 
Research at MSRI is supported in part by NSF grant DMS-9701755. 



References 

[1] M. Mehta, Random matrices, Academic Press, New York, 1991. 

[2] D.V. Voiculescu, K.J. Dykema, A. Nika, Free random variables, CRM Monogr.Ser.l, AMS., 

Providence, 1992. 
[3] B. Khoruzhcnko, L. Pastur, A. Stojanovich, (in preparation). 

[4] N. Akhiezer, I. Glazman, Theory of linear operators in Hilbert space, Frederick Ungar, New 
York, 1961. 

[5] D. Ruelle, Statistical mechanics, Benjamin, New York, 1969. 

[6] L. Pastur, A. Khorunzhy, On the eigenvalue distribution of random matrices with independent 
entries, Preprint N21 of the Institute of Low Temperature Physics, Kharkov, 1989. 

[7] A. Khorunzhy, S. Molchanov, L. Pastur, On the eigenvalue distribution of band random matri- 
ces, Teor. Math. Phys., 90, 108-118 (1992). 

[8] A. Khorunzhy, B. Khoruzhcnko, L. Pastur, M.Shcherbina, Large-n limit in the statistical me- 
chanics and the spectral theory of disordered systems, In: Phase Transitions and Critical Phe- 
nomena,]^ C.Domb, J.Lebowitz (Eds.), Academic Press, New York, 1992. 

[9] A. Khorunzhy, L. Pastur, Limits of infinite interaction radius, dimensionality and the number 
of components for random operators with off-diagonal randomness, Commun. Math. Phys., 
153, 605-646, (1993). 

[10] B. Khoruzhenko, W. Kirsch, L. Pastur, The interband light absorption coefficient in the weak 
disorder regime: an asymptotically exactly solvable model , J. Phys. A: Math. Gen., 27, 2527- 
2543 (1994). 

[11] A. Khorunzhy, Eigenvalue distribution of large random matrices with correlated entries, 

Mathem.Phys., Analysis and Geometry, 3, 80-101 (1996). 
[12] J.-L. Pichard, Random transfer matrix theory and conductance fluctuations , In: Quantum 

Coherence in Mesoscopic Systems, B.Kramer (Ed.), Plenum, New York, 1991. 
[13] R. Muirhead, Aspects of multivariate statistical theory, Wiley, New York, 1982. 
[14] V. Marchenko, L. Pastur The eigenvalue distribution in some ensembles of random matrices 

Math.USSR Sbornik, 1, 457-483 (1967). 
[15] D. Bessis, C. Itzykson, J. Zuber, Quantum field theory techniques in graphical enumeration, 

Adv.Appl.Math., 1, 109-157 (1980). 
[16] F. Xu A random matrix model from two-dimensiona Yang-Mills theory, Commun. Math. Phys., 

190, 287-307 (1997). 
[17] D. Voiculescu (Ed.) Free probability theory, AMS, Providence, 1997. 



24 



LEONID PASTUR 



[18] P. DiFrancesco, P.Ginsparg, J. Zinn- Justin, 2D gravity and random matrices, Phys. Reports, 
254, 1-133 (1995). 

[19] A. Boutet de Monvel, L.Pastur, M. Shcherbina, On the statistical mechanics approach to the 
random matrix theory: the integrated density of states, J. Stat. Phys., 79, 585-611 (1995). 

[20] L.Pastur, M. Shcherbina, Universality of the local eigenvalue statistics for a class of unitary 
invariant ensembles, J. Stat. Phys., 86, 109-147 (1997). 

[21] N.I. Muskhelishvili, Singular integral equations, Nordhoff, Groningcn, 1953. 

[22] L. Pastur, On the universality of the level spacing distribution for some ensembles of random 
matrices, Lctt.in Math.Phys., 25, 259-263 (1992). 

[23] L. Pastur, V. Vasilchuk, On the law of addition of random matrices, (in preparation) 

[24] A. Zee, Law of addition in random matrix theory, Nucl.Phys., B474, 726-744 (1996). 

[25] R. Speicher, Free convolution and the random sum of matrices, Publ.RIMS, 29, 731-744 (1993). 

[26] P.Neu, R. Speicher, Rigorous mean-field model for coherent-potential approximation: Anderson 
model with free variables, J.Stat.Phys., 80, 1279-1308 (1995). 

[27] V. Vasilchuk, On the law of multiplication of random matrices, (in preparation) 

[28] L. Pastur, On the spectrum of random matrices, Teor. Math. Phys., 10, 67-74 (1972). 

[29] V. Girko, Theory of random matrices, Kiev, Vyshcha Shkola, 1975 (in Russian). 

[30] V. Girko, Theory of random matrices Kluwer Academic Publishers (1990). 

[31] A. Khorunzhy, B. Khoruzhenko, L. Pastur, Random matrices with independent entries: asymp- 
totic properties of the Green function, J. Math. Phys., 37, 5033-5060 (1996). 

[32] L. Pastur, Eigenvalue distribution of random matrices, Ann.l'Inst.H.Poincare,64, 325- 
337(1996). 

Department of Mathematics, University Paris-VII, 2, place Jussieu, 75251, Paris, 
France 

E-mail address: pastur@math.jussieu.fr 

Mathematical Division, Institute for Low Temperatures, 47, Lenin's Ave., 310164, 
Kharkov, Ukraine 

E-mail address: lpastur@ilt.kharkov.ua 



